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ABSTRACT 


The  Buffon  needle  probleoi  and  soae  variations 
are  used  to  Illustrate  classical  statistical  Methods 
of  estimation  and  to  lead  Into,  and  contrast  with, 
the  problems  which  arise  when  a sample  of  some  random 
structure  Is  the  data.  The  flavor  of  these  problems 
Is  conveyed  largely  by  discussion  of  the  simplest, 
and  most  described,  case,  that  of  point  processes. 
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1.  INTRODUCTION 


Th«  Buffon  needit  problen  h«s  Its  origins  In  gambling.  But,  un- 
like card  and  dice  games.  It  requires  some  geometry  and  calculus. 

Thus  It  was  a bold  step  In  generalizing  the  Idea  of  mathematical 
probability.  The  needle  problem  and  minor  variants  have  been  studied 
by  many  authors  up  to  the  present  day  as  a means  of  estimating  v 
statistically.  More  sensibly,  however,  they  provide  a means  of  In- 
ferring the  ratio  of  the  length  of  the  needle  to  the  scale  of  the 
regular  network  onto  which  It  Is  thrown.  And  of  course  this  Is  the 
origin  of  modern  methods  for  the  sampling  study  of  geometric  bodies. 
Thus  It  seems  most  appropriate  at  this  Symposium  to  begin  my  talk  by 
going  over  this  work  which  used  only  quite  standard  statistical  Ideas. 

However  when  the  needle  Is  replaced  by  a more  general  probe  or 
sampling  window  --  In  other  words,  some  set  — and  the  regular  net- 
work Is  replaced  by  a random  structure,  we  must  face  statistical  prob- 
lems of  quite  a different  character.  They  are  akin  to  those  met  In 
time  series  analysis  which  has  a long  history  by  statistical  standards 
but  In  which  there  Is  still  some  confusion  between  the  exploratory, 
modelling  and  confirmatory  aspects  and  a lack  of  communication  between 
probablllsts  and  statisticians.  The  spatial  problems  are  decidedly 
more  difficult.  There  are  not  so  many  explicit  probability  models  to 
get  statistical  experience  with.  One  must  consider  the  shape  as  well 
as  the  size  of  sampling  windows.  The  data  has  a more  awkward  form. 

The  second  part  of  the  lecture  will  therefore  merely  give  some  Idea 
of  these  problems  and  their  literature.  Hopefully  other  speakers 
will  address  them  In  detail.  But  It  Is  clear  that  statistical  geo- 
metry or  morphology  Is  Just  beginning. 


2. 


I 


1 1 

i*< 


l.  BUFfOW  HOBIEMS 

The  clessicti  problen  considers  a parallel  grid  with  spacing 
and  a needle  of  length  t where  £ < a.  If  the  needle  Is  tossed 
that  Its  position  and  orientation  are  random,  the  probability  p 
It  cuts  a grid  line  Is  given  by 

w a, 

Uspensky  (1937).  Kendall  and  Moran  (1963),  for  example,  give  proofs  of 
this  and  most  of  the  results  used  below. 

Define  d * r • <*  * f*  suppose  that  n Independent  trials 
(tosses)  yield  C cuts.  Then  jC(C)  ■ Binomial  (n,p)  where  p ■ 2^r. 
The  statistical  problems  that  arise  are 

(1)  estimate  p,  I.e.,  estimate  v 

(11)  estimate  r,  I.e.,  estimate  a If  £ Is  known 

(111)  test  that  r • r^  (known). 

The  most  complete  reference  on  (1)  In  this  case  and  those  given  below 
Is  Perlman  and  Michura  (1975).  Oddly  no  one  seems  to  have  considered 
(11)  and  (111).  Trivial  though  they  are,  they  are  the  prototypes  of 
the  real  problems. 

The  likelihood  of  the  data  when  c cuts  are  observed  Is 

(")P^(1-P)''~^  snd  the  number  of  cuts  Is  a complete  sufficient  statis- 
tic for  p.  Thus  we  may  assert:  among  all  functions  f(c)  such  that 
Ef(C)  ■ p,  for  all  p In  (0,1),  It  Is  true  that 

var  f(C)  I var  £ - n p(l-p)  . (2) 

Thus  the  obvious  estimator  of  p,  p ■ C/n  Is,  uniformly  In  p,  the 


a 

so 

that 

(1) 


I 
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minimum  variance  unbiased  estimator  of  p. 

If  one  knows  r and  wishes  to  estimate  v,  the  same  Is  true  of 


However  - l/$  Is  biased.  If  n Is  large,  we  may  argue  that 


w as  n . 


Further  It  follows  from  (3)  that 


Thus  If  one  Is  set  on  estimating  ir  this  way  and  can  "design  the 
experiment”  one  should  take  t * a,  since  this  choice  minimizes  (4). 
In  this  case  var  = 5.63/n.  Lazzerinl  (1901)  conducted  such  an 
experiment  with  n ■ 3408  and  found  - ir  ■ 3 » 10*^.  As  Kendall 
and  Moran  (Ibid.)  suggest,  he  must  have  stopped  when  he  noticed  the 
remarkable  and  fortuitous  accuracy! 


In  fact  we  know  4 and  are  more  likely  to  want  to  know  a.  By 

the  same  arguments  used  above,  a • , and  we  will  have 

c 

Ea  a as  n • and 


4. 


Again  If  Me  can  choose  the  needle  size  we  should  try  to  nake  It 
near  to  but  less  than  a.  Thus  the  effort  to  optimize  may  lead  to 
bias  since  (1)  Is  false  when  t > a.  Trivially  If  n Is  large,  a 
Is  Gaussian  so  tests  are  easy  to  make. 

Instead  of  a parallel  grid,  Laplace  considered  a rectangular 
grid,  the  A lines  being  a apart,  the  B lines  being  b apart. 
He  showed  that  the  probability  that  the  needle  cuts  at  least  one  line 


2^(a-»-b)  - V 
irab 


a fascinating  formula  whose  direct  derivation  Is  tricky  so  that  It 
seems  easier  to  get  It  from  a more  general  Crofton  argument.  To 
Illustrate  my  points  here  set  a.  * b and 


^ ^ - r 

a “ F * " • 


Introducing  the  notation 

■ Prob(need1e  cuts  a B line  but  not  an  A line), 

Pj^  • Prob(need1e  cuts  neither  an  A nor  a B line), 

etc.,  formula  (6)  Is  clearly  1 - P^  . Further  If 

P^  ■ Prob(need1e  cuts  an  A line) 

Pg  • Prob(need1e  cuts  a B line) 


then 


'a  ■ 'ab  * • 

'b  • 'aB  * ’’iB  • ’ 

'ab  * ’’aF  * 'Kb  * ’’ff  * *• 


5 


These  three  equations  plus  (6)  for  1 - yield 

'aB  ■ • ''SB  ' ''aF  ’ <') 

and  for  brevity  we  follow  Perlman  and  Ulchura  In  writing 

Pjg-  - 1 - (4r-r^)^  » 1 - m*  . (8) 


Let  n trials  yield  results  (In  an  obvious  notation) 

N - (n^g,  n^y,  n^g.  nj^} . Then 

J^(N)  • 4->nom1al  (n;  P^g,  P^g-,  Pjjg,  P^)  . 

Thus  the  likelihood  of  the  data  Is  proportional  to 

"AB  "aF  "JiB 

^ • '’ab  '’JTb 

Defining 

Nq  * nj^  ■ I no  cuts 

"l  “ "aB-  * ^ 

N,  ■ n^g  ■ # 2 cuts 

n « Nq  ♦ Nj  ♦ N2 

and  using  (7)  and  (8),  L may  be  written  as 


Nn  N,+N,  2N,+N,  N, 

(l-m4)  ° 1 2 ^ 2 1 1 


(9) 


Thus  If  r Is  known.  Nq  or  N2  Is  a complete  sufficient 

statistic  for  i and  i(N^-i-N2)  ■ B1nom1a1(n»  m4)  so  the  story  of  the 
estimator  of  n follows  the  previous  pattern.  The  resulting  est1ma> 
tor  *2  1*rge»  a variance  equal  to  0.47/n  so  that  the 

extra  work  In  using  a square  grid  yields  an  estimator  which  Is  12 
times  as  efficient  as  that  for  the  parallel  grid. 


6. 


HoHtver  to  estimate  r knowing  ir  Is  quite  different.  The 
practical  method  Is  to  choose  r to  maximize  the  likelihood  (9). 
Setting  3 log  L/3r  equal  to  zero  leads  to  the  equation 


2r-4 


ZMg+Nj 


1 - 


4r-r* 


HL. 

2-r 


(10) 


which  must  be  solved  Iteratively  to  yield  r.  The  standard  theory  of 
maximum  likelihood  estimation  gives  us  an  asymptotic  formula  for  r. 


var  f - 


log  t 


3r 


r 


and  asymptotic  normality  of  f so  that  tests  can  be  made. 
The  rectangular  grid  follows  the  same  pattern  with 


(ii: 


" a • *^2  * F * *■  * ^“ootlon  of  r^  and  r^  which  are  esti- 

mated by  solving  3 log  L/3rj  ■0.  3 log  L/3r2  * 0.  Other  regular 
networks  do  not  Introduce  the  need  for  further  techniques. 

Above  we  considered  only  the  case  of  £ $ a « b.  The  case  when 
t Is  much  greater  than  a Is  simple  and  Instructive  to  consider. 
Let  then  this  long  needle  Intersect  the  B lines  at  an  angle  0. 
Then  6 Is  uniformly  distributed  on  (0,ir/2).  If  we  define 

• I A lines  cut  s r cos  6 • 

Ng  ■ I B lines  cut  s r sin  8 • 

N ■ # lines  cut  : r(cos  0 * sin  0)  » 

we  have 


7. 


r cos  e de  - 


^ - £(«,)  . 


E(r^  > 2r^  cos  9 sin  0)  , 

r^d  ♦ |)  . 


va  r fl  * r 1 + 


j-©')  ■ 


Mence  If  we  make  n throws  and  find  Z as  the  average  number  of 
ruts.  It  will  be  an  unbiased  estimator  of  E(N).  Hence  the  cstlmatcr 
of  V that  Is  suggested,  following  our  ealler  work,  is 


3 — 

^ C 


and  an  easy  calculation  shows  that 


var 

3 16r^  " 


. 0.C095 


While  this  seems  a great  Improvement,  we  will  show  below  that  one  can 
do  better  still  with  this  experiment. 

To  estimate  r from  c , the  "natural*  method  Is  to  set  r ■ 


* ,,2  2 f M 

with  var  estimate  too  cin  be  Improved 

because  neither  Is  the  maximum  likelihood  estimator,  as  was  true  of 
our  first  three  examples. 


•Hw 
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From  (12),  N = /?  r cos  le  - so  that 


Prob(N  < k)  * Prob(cos  < k/r/T)  where  Is  uniform  on 

(-it/4,  it/4).  Let  ■ COS*'  . Then 

° r/? 


Prob(N  < k)  * 2 Prob(i|»  > *q)  , 

(13) 


so  that  the  probability  density  of  N at 
tive  of  (13)  with  respect  to  k,  namely 


k 


Is  the  partial  deriva- 


Thus  given  counts  k, k in  n trials,  their  likelihood  is 

in 


provided  r s all  the  k^'s  < /?  r,  and  zero  otherwise.  Thus  the 
maximum  likelihood  estimate  of  r is 


r*  ■ ^ max(kj, . . . ,k^)  , (14) 

not  r ■ irc/4.  It  may  be  shown  that  the  variance  of  r*  is  of  order 
n"^,  not  order  n"'  like  that  for  r.  Thus  for  large  n,  r*  is  a 
very  much  better  estimator  than  r.  This  shows  dramatically  that  the 
usual  practice  in  geometrical  statistics  of  obtaining  estimators  by 
equating  theoretical  and  observed  means  may  be  very  inefficient.  So 
much  for  the  "long  needle."  Other  details  may  be  found  in  Diaconis 
(1976)  and  in  a forthcoming  monograph  by  H.  Solomon. 
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Buffon's  needle  may  be  used  to  obtain  a connection  with  a quite 
different  aspect  of  geometrical  statistics.  Let  us  analyze  the  toss- 
ing of  the  needle.  Suppose  now  that  an  origin  Is  marked  on  one  of  th 
lines  of  the  parallel  grid  and  that  the  needle  Is  thrown  so  that  Its 
center  rests  on  the  plane  a distance  X from  the  marked  line.  Let 
^(X)  « Gaussian  (O.o^)  so  It  has  probability  density 

f(x)  • {a/2¥) "^exp( -x^/2a^)  . 

It  is  then  clear  that  If  Y Is  the  distance  from  the  center  of  the 
needle  to  the  nearest  line  below  It, 


Prob{x  < Y s X + dx)  - J f(x  - Vtt)dx 

■ g(x}dx  , say. 


(15) 


The  density  g(x)  Is  concentrated  on  (0,a)  and  g(x)  Is  periodic, 
period  a.  Thus  we  may  write 


S(x)  " I g^exp(-2TTl jx/a)  . (16) 

It  is  shown  that  (see  Hartman  and  Watson  (1974)]  this  density  can  be 
very  well  approximated  by 

(2iTlQ(tc))  “'exp  te  cos  (2irx/a)  (17) 

where  k is  a suitably  chosen  function  of  a and  that  by  certain 
randomizing  an  exact  result  may  be  obtained.  One  of  the  commonest 
distributions  for  describing  non-uni formly  distributed  angles  (which 
we  would  need  If  we  wished  to  give  the  needle  a preferential  orienta- 
tion) Is  the  von  Mises  distribution 


I 

I 


10. 


(2fflQ(t<:))  “^exp  K cos  0 


(18) 


Our  final  two  examples  lead  into  random  structures.  Suppose  tha 
the  spacings  of  the  parallel  grid  are  Identically  and  Independently 
distributed  (I.I.D.)  with  some  density  function  h{a)  wiilch  Is  zero 
when  CL  < a^.  Consider  a needle  of  length  £ < aQ  tossed  at  random. 
Then 


Prob(center  of  the  needle 
falls  In  a space, 
a < space  < a + da)  » 


ah(a) da 


C 


ah(a) da 


Since 


ah(a)da 

E(a) 


Prob(cut| space  a) 
Prob(cut) 


2 I 


fl  L . 

J ff  a E(a} 
‘>0 

2 I 
^ * 


Thus  by  repeating  this  experiment  the  only  thing  we  can  learn  about 
the  spacings  is  E(a).  There  Is,  for  example,  no  way  one  can  check 
whether  they  are  I.I.D.  This  would  require  a long  needle. 

Thus  let  us  consider  an  Infinitely  long  needle  and  suppose  that 
we  could  know,  after  It  Is  tossed  at  random  onto  an  arbitrary  para- 
llel grid,  the  sequence  of  spaces  on  the  needle  between  line  cross- 
ings, {s^}.  If  the  needle  makes  an  angle  e (which  we  do  not  know) 
and  the  grid  spacings  are  then  • s^  sin  6,  for  all  positive 


11. 


and  negutive  Integers  1.  From  only  the  sequence  (s^)  we  can  check 
ail  properties  of  the  (a^}  sequence  that  do  not  depend  upon  scale, 
e.g.,  that  It  Is  stationary,  etc. 

The  cut  points  on  the  needle  form  a Point  Process  In  one  dlsisn- 
slon.  The  discussion  of  Point  processes  In  space  Is  the  largest 
aspect  of  geometric  probability  and  statistics  In  the  modem  sense, 
the  topic  to  which  we  now  turn. 

3.  WHAT  IS  STATISTICAL  GEOMETRY? 

Everyone  is  fairly  clear  what  Is  meant  by  geometry  but 
statistics  Is  less  well  defined.  It  has  a number  of  facets  • 

exploring  dcta  for  regularities,  I.e.,  patterns;  (11)  estimating 
"population"  characteristics  from  a sample;  (111)  testing  hypotheses: 
(1v)  designing  sampling  plans  to  be  effective  and  efficient  (usually 
by  Including  a random  element). 

Probability  models  enter  (1)  to  (1v)  In  several  ways: 

a)  by  a scientific  mechanism  or  model, 

b)  by  assumption, 

c)  via  a random  sampling  plan, 

and  to  different  degrees.  In  (1)  they  nay  not  enter  explicitly  at 
all. 

In  statistical  geometry  our  data  will  be  a sample  from  some  geo- 
metrical "population.” 

Such  definitions  do  not  convey  much,  so  we  now  give  some  examplr 
of  probloas  and  the  groups  that  pursue  them.  (A)  Grenander's  books 
on  "Pattern  Synthesis*  (1976)  cover  a vast  area  In  a novel  way  not 
represented  at  all  at  this  Symposium,  and  1 think  they  are  of  basic 
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Importance.  He  has  developed  an  abstract  way  of  generating  and  distort- 
ing patterns  and  then  restoring  them.  The  latter  Is  of  course  statis- 
tical. He  gives  a wealth  of  diverse  examples;  one  of  the  simplest 
Is  discussed  below.  (B)  Classical  problems  such  as  may  be  found  In 
the  Kendall  and  Horan  (1963)  book.  (C)  The  publications  of  the 
Fontainebleau  School  of  Mathematical  Morphology  represent  a different 
line  again.  Their  major  achievement  seems  to  me  to  be  the  wedding  of 
the  Image  analyzer  and  mathematical  description  of  the  objects 
scanned.  While  much  practical  work  Is  done,  the  publications  deal 
more  with  the  mathematical  theory  than  with  the  statistical  aspects. 

(D)  The  Stochastic  Geometry  pursued  In  Cambridge  by  0.  6.  Kendall 
(see,  e.g.,  Harding  and  Kendall,  1974)  and  others  overlaps  theoreti- 
cally with  the  French  School  but  has.  It  seems,  purely  mathematical 
motivations.  Like  so  much  of  this  literature.  It  has  not  been  re- 
duced to  a level  of  mathematical  simplicity  for  practical  statistical 
use.  (E)  The  Point  Process  literature,  stemming  from  Bartlett  (see, 
e.g.,  1963,  1964,  1976)  originated  In  practical  statistical  problems 
and  mainly  In  one  dimension.  It  is  now  pursued  at  a highly  mathe- 
matical level  by  Krickeberg  (1977)  and  other  Europeans  In  many  dimen- 
sions. Earlier  practical  work  In  Forestry,  especially  Natem's  (2960) 
has  led  to  many  papers  — see,  e.g.,  the  Issuesof  Blometrlka. 

Ripley's  recent  papers  (1977a,  1977b)  have  a combination  of  theory 
and  practice  and  extensive  bibliographies. 

While  there  are  many  mechanisms  for  generating  point  processes 
In  time,  the  few  that  do  so  In  space  are  summarized  by  Ripley  (1977a), 
for  example.  The  main  emphasis.  In  line  with  second  order  stationary 
processes.  Is  the  definition  and  estimation  of  functions  that  control 

• I 
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the  enhancement  or  Inhibition  of  neighboring  points.  Here,  as  In  time 
series  analysis*  there  Is  a large  exploratory  element.  Even  If  there 
were  parametric  models*  It  would  rarely  be  possible  to  write  down  the 
likelihood  of  the  data  so  that  the  t1me*honored  statistical  methods 
Illustrated  earlier  cannot  be  used.  Computers  are  essential  for  al- 
most all  calculations*  unlike  the  Buffon  problems*  e.g.«  variances 
must  usually  be  found  by  simulating. 

In  practice  we  will  often  want  random  sets*  rather  then  the 
Poisson  fields  of  points*  lines*  flats*  etc.  that  are  most  often 
discussed.  In  his  1967  book  Matheron  made  one  of  the  early  models 
chat  can  be  dealt  with  easily  — the  Boolean  scheme.  Here  I.I.O. 
copies  of  a random  set  are  attached  to  a Poisson  field  of  points 

{x^jx^eX}  In  a vector  space  to  obtain  the  random  set 

A ■ U *<)•  Such  a set  Is  Intuitively  stationary*  1.e.» 

x^X  ’ ’ 

spatially  homogeneous  though  not  necessarily  Isotropic.  If  we  know* 
for  any  fixed  set  B* 

q(B)  • ProbCkOB  ■ 0} 


then 

ProbCADt  • ()  • e>p(-if  (1-4(B*0)AU  (») 

V 

where  X Is  the  Intensity  of  the  Poisson  process.  In  this  Symposium 
Coleman  went  further  In  this  construction  than  (19)  which  Is  the  zero 
term  of  a Poisson  distribution. 

Time  series  analysis  Is  about  100  years  old.  It  began  at  a 
practical  endeavor*  became  very  mathematical  and  It  Is  only  recently 
that  practical  books  and  programs  have  been  readily  available.  The 
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time  lag  for  this  subject  could  be  greatly  shortened  If  theoreticians 
would  make  the  effort  to  write  for  practical  users  and  not  only  for 
other  mathematicians. 

We  conclude  with  the  simplest  Instance  of  problems  In  Grenander's 
book  (Ibid.).  It  Illustrates  (1)  how  a finite  window  Is  different 
from  a finite  sample  of  I.I.O.  observations*  (11)  the  use  of  Fourier 
analysis.  It  Is  the  restoration  of  a linear  lattice  whose  points 
have  been  Independently  displaced.  The  complete  set  of  points  Is 

■ a ♦ v5  ♦ n^  (v  ■ ....  -1,  0,  1*  ...) 

where  a Is  a phase,  ( ■ the  lattice  spacing  (unknown)  and  n^  Is 
the  noise.  The  window  Is  the  Interval  (O.L).  When  the  noise  Is  small 
with  respect  to  L»  almost  all  the  points  that  should  be  In  (O.L)  will 
be  there  and  no  two  points  will  have  their  true  order  Inverted.  Then 
we  have  an  ordinary  regression  problem  In  estimating  a and  C. 

When  the  noise  Is  not  small,  successive  X points  may  not  have  suc> 
cessive  Indices  and  the  "wrong"  points  may  be  In  the  window  -•  this 
Illustrates  point  (1).  (This  model  Is  essentially  the  same  as  that 
set  up  by  0.  G.  Kendall  (1974)  to  detect  a unit  of  measurement  In  an 
archeological  site.) 

Here  one  automatically  thinks  of  Fourier  analysis.  To  save  time, 
set  a ■ 0 and  define. 

♦(w)  • f I exp  1 uX  . (20) 

^ X^e(O.L)  '' 

m(w)  • E4(u)  . 

1 

■ ^ J^exp(1ux)p(x)dx. 


j 

I 
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■ I f (x  - v?) , 

V 

• density  of  the  noise  n^ 

• Z Pi,  exp( -27rik/C) . 

expQw  - - I 

- Z p.  . 

" iiu  - 

0,  u f multiple  of  2vk/C 
p|^.  u ■ multiple  of  2«k/C  • 

Thus  we  would  hope  to  see  « pattern  of  peaks  near  the  points  2irk/C 
from  which  we  would  first  try  to  see  i_f  there  is,  a pattern*  and  If  so 
to  estimate 

var  - r 1 

where  f*  Is  tho  Fourier  transform  of  f so  that 


Thus  var(f(w))  may  also  help  us  learn  about  the  noise  since  this  Is 
governed  by  f. 

The  same  computation  with  a Poisson  process  leads  to 

C(d)  ■ X ^ t • very  different  picture  also 

seen  with  all  renewal  processes.  Thus  # does  not  differentiate  be- 
tween stationary  point  processes  but  one  hopes  that  the  variance  migh 
If  we  define 


where 

p(x) 

f(x) 

and 

p{x) 

Then 

m((i>) 

As  L 

m(fc)) 
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Me  have 


N(x)  * I points  In  (0.x)  , 
♦ (w)  ■ |^exp(1ux)dN(x)  , 


Assune  Mith  Bartlett  (1976)  that 

,2 


E(dH(x))‘  - Xdx 
EdN(x)dN(y)  - (X^  * H(x-y)}dxdy  . 


(21) 


If  U(*)  In  (21)  Is  Identically  zero,  the  points  are  Poisson.  Let 
U(v)  - W(-v).  It  Is  clear  that  If  W(v)  Is  positive,  a point  at  y 
«eans  that  there  1$.  relative  to  the  Poisson  process,  acre  chance  of 
having  a point  at  y v.  I.e..  enhanceeent.  Negative  values  aean 
Inhibition.  Now 


(22) 


which  verifies  our  notion  that  knowledge  of  |4|‘  should  yield  In- 
formation about  the  function  W(v). 

In  this  use  of  Fourier  analysis  one  should  note  that  the  F.F.T. 
cannot  be  used;  It  Is  hard  to  adjust  for  bias  and  finite  L and  hard 
to  find  the  variance  of  |d)^.  even  In  R'.  The  vagueness  In  these 
last  two  paragraphs  Is  to  some  extent  unavoidable.  When  dealing  with 
unknown  functions,  one  simply  has  to  use  Judgment,  try  various  tricks 
with  the  computer  --  there  cannot  be  any  simple  and  apparently  clear 
cut  methods  of  the  t-test  type. 


17. 
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some  random  structure  Is  the  data.  The  flavor  of  these  problems 
Is  conveyed  largely  by  discussion  of  the  simplest,  and  most  des- 
cribed, case,  that  of  point  processes.^ 
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